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ABSTRACT. It is shown that the hyperbolic functions can be associated 
with selfdecomposable distributions (in short: SD probability distributions 
or Levy class L probability laws). Consequently, they admit associated back- 
ground driving Levy processes Y (BDLP Y). We interpret the distributions 
of y(l) via Bessel squared processes, Bessel bridges and local times. 
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1. Introduction and terminology. The aim of this note is to pro- 
vide a new way of looking at the hyperbolic functions: cosh, sinh and tanh, 
or their modifications, as the members of the class SD of selfdecomposable 
characteristic functions (often called class L, after Paul Levy). Analytically 
one says that a characteristic function <fi is selfdecomposable, we simply write 

<p e SDjf 

V(0<c< 1)3 p c V(t€K) <j>{t) = <p(ct) Pc (t), (1) 

where p c is also a characteristic function. Let us recall here that class SD 
is a proper subset of ID, the class of all infinitely divisible characteristic 
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functions, and that the factors p c in (1) are in ID as well ; cf. Jurek and 
Mason (1993), Section 3.9., or Loeve (1963), Section 23 (there this class is 
denoted by < Jl). We also will use the convention that a random variable X 
(in short: r.v. X) or its probability distribution fix or its probability density 
fx is selfdecomposable if the corresponding characteristic function is in the 
class SD. Furthermore, the equation (1) describing the selfdecomposability 
property, in terms of a r.v. X means that 

X G SD iff V(0 < c < l)3(r.v. X c ) X = cX + X c , 

where the r.v. X and X c are independent and = means equality in distribu- 
tion. 

For further references let us recall the main properties of the selfdecom- 
posable distributions (or characteristic functions or r.v's): 

(a) SD with the convolution and the weak convergence forms a closed con- 

volution subsemigroup of ID; 

(b) SD is closed under affine mappings, i.e., for all reals a and b one has: 

G SD iff e ibt (j)(at) G SD. 

(c) X G SD iff there exists a (unique) Levy process Y(-) such that 

X = J oo e-W(s), wh ere Y is called the BDLP (background driving 
Levy process) of the X. Moreover, one has that E[log(l + |y(l)|)] < oo. 

[ID\ og will stand for the class of all infinitely divisible laws with finite 
logarithmic moments.] 

(d) Let (f) and ip denote the characteristic function of X and Y{1), respec- 

tively, in part (c). Then one has 

log 0(0 = jjlog^f , i.e., m = exp[t(lo g 0(t))'],t ^ 0,#)) = 1. 

(e) Let M be the Levy spectral measure in the Levy-Khintchine formula of 

4> G SD. Then M has a density h(x) such that xh(x) is non-increasing 
on the positive and negative half-lines. 

Furthermore, if h is differentiable almost everywhere then dN(x) = 
— (xh(x) dx is the Levy spectral measure of ip in (d). 
Finally, one has also the following logarithmic moment condition 
If\x\>e} 1°S(1 + \x\)dN(x) < oo, for all positive e. 
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Parts (a) and (b) follow directly from (1). For (c) and (d) cf. Jurek & 
Mason (1993), Theorem 3.6.8 and Remark 3.6.9(4). Part (e) is Corollary 1.1 
from Jurek (1997). 

In this note we will characterize the BDLP's (or the characteristic func- 
tions ip in d) for the hyperbolic characteristic functions. The main result 
shows how to interpret these distributions in terms of squared Bessel bridges 
(Corollary 2) and squaredd Bessel processes (Corollary 3). 

2. Selfdecomposability of the hyperbolic characteristic func- 
tions. For this presentation the most crucial example of SD r.v. is that 
of the Laplace (or double exponential) random variable rj. So, rj has the 
probability density \e~\ x \x E R, and its characteristic function is equal 



1 r r°° p~\ x \ 

« i) = _ = exp[y (e --l) 



oc 
OO 



dx 

X 



r* r r°° i dv 

= exp / / [e lvx - l)e-^dx — e SD. (2) 
Jo l J-oo \ v 

To see its selfdecomposability property simply note that 
<pJt) 1 + cH 2 ~ .„ 1 . , , . . , 

—7 — - = 7T = c 1+ 1— c ^ is the characteristic function otp c , 

(f) v (ct) 1 + t 2 y J l + t 2 1 ' 

in the formula (1). The rest follows from appropriate integrations; cf. Jurek 
(1996). 

Another, more " stochastic" argument for selfdecomposability of Laplace 
rv rj, as a counterpart to the above analytic one, is as follows. 

Firstly, notice that for three independent rvs S(l), S(l) and b c , where the 
first two have exponential distribution with parameter 1 and the third one 
has Bernoulli distribution (P(b c = 1) = 1 — c and P{b c = 0) — c), one has 
equality 

£{l) = c£{l) + b c £{\), 

which means that £(1) is a selfdecomposable rv. (The above distributional 
equality is easily checked by using the Laplace or Fourier transform.) 

Secondly, taking two independent Brownian motions B t , B t , t > and 
independently of them an exponential rv £(1) satisfying the above decompo- 
sition, we infer that 

BeiX) = VcB £{ i) + B b j {1) , 
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and thus proving that stopped Brownian motion -B^(i) is selfdecomposable 
as well. 

Thirdly, let us note that -B^(i) has the double exponential distribution. 
More explicitly we have 

1 ~~ \~ t 

For more details and a generalization of this approach cf. Jurek (2001), 
Proposition 1 and Bondesson (1992), p. 19. 

PROPOSITION 1. The following three hyperbolic functions: — ^r, "77^77; 
ta ° h * , t G R, are characteristic functions of selfdecomposable probability dis- 
tributions, i.e., they are in the class SD. 



Proof. From the following product representations : 

-IV^ )' sinhz = zj](^ 

fc=l v ' k=l 



cosh^ = J](l+ ), sinhz = zj](l + -^), (3) 

i.— 1 \ / i— i 



for all complex z, and from (2) with (a) we conclude that the first two 
hyperbolic functions are characteristic functions from SD. Moreover, these 
are characteristic functions of series of independent Laplace r.v.; cf. Jurek 
(1996). 

Note that for < a < b the fraction 

1 _|_ a 2 t 2 a 2 a 2 l 

7777 = 77 1 + (1 - 77 h 7777 is a characteristic function, 

1 + b z V tr cr 1 + 01 

tanht fr l + (kny 2 t 2 



and so is 



t _ n 1+((A; _i )7r) -2 t2 ' 



as a converging infinite series of characteristic functions of the above form. Its 
selfdecomposbility follows from Yor (1997), p. 133, or Jurek (2001), Example 
1(b). □ 

REMARK 1. The selfdecomposability of ^f 1 , i.e., the formula (1), would 
follow in an elementary manner if for all < c < 1 and all < w < u, the 
functions 

1 + (c 2 u + w)t 2 + c 2 uwt 4 1 + wt 2 1 + c 2 ut 2 

; 5 — — — = — - were characteristic functions. 

1 + (u + c 2 w)t 2 + c 2 uwt A 1 + ut 2 1 + c 2 wt 2 
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[The above is a ratio of two fractions of the form as in the product repre- 
sentation of tanht/t, with the fraction in the denominator computed at ct.\ 
However, they can not be characteristic functions ! Affirmative answer would 
mean that Laplace rv is in L\ (these are those SD rv for which BDLP Y{1) 
in (c), is SD). Equivalently, the characteristic function p c in (1) is in SD).) 
But from (2) we see that Y(l) for Laplace rv r] has compound Poisson distri- 
bution with Levy spectral measure dM(x) = e~^dx which does not satisfy 
the criterium (e). Cf. also Jurek ( 1997). 

3. The BDLP's of the hyperbolic characteristic functions. Since 
the three hyperbolic characteristic functions are infinitely divisible one can 
insert them into Levy processes: C s , S s , T s , for s > 0, corresponding to cosh, 
sink and tank characteristic functions. Those processes were studied from 
the ID class point of view in the recent paper Pitman- Yor (2003). Here we are 
looking at them from the SD class point of view, i.e., via the corresponding 
BDLP's. 

Below (j) with subscript C, S or T denotes one of the three hyperbolic 
characteristic functions, M with similar subscripts denotes the Levy spectral 
measure in the appropriate Levy-Khintchine formula, furthermore ifj with the 
above subscripts is the corresponding characteristic function in the random 
integral representation (properties (c) and (d) of class SD) and finally N 
with one of the above subscripts is the Levy spectral measure of ip (as in 
(e)). Thus we have the equalities : 

Mt) = Mt)-M), i-e., ^ = -^-^ (4) 



M d {-) = M § (-) + M t (•), where dM c^ _ 1 



dx 2x sinh(7ra;/2) ' 

dMa(x) e -^\n \ -\ x \ 

— = n t = -TT cot h (— ) - 1); 

dx 2a;sinh(7ra;/2 2\x\ K y 2 ' h 

dM f (x) 1 e-^l^l/ 4 1 M , . . . /r , 

1^2- = —— — = -— i - tanh tt x A )]. 5 

dx 2\x\ cosh(vr|x|/4) 2\x\ 1 y ' 

These are consequences of the appropriate Levy-Khintchine formulas for the 
hyperbolic characteristic functions or see Jurek (1996) or Pitman- Yor (2003) 
or use (3) and the product formulas for cosh z, sinh z; (for tanh z use the ratio 
of the two previous formulas). 
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COROLLARY 1. For the SD hyperbolic characteristic functions 0<a 7 0^ 
and 4>f, their background driving characteristic functions are ip^, ipg and ipf, 
where 



^c(^) = ^s^'^fit)' = exp[— t tanht], tp§(t) = exp[l— t cotht], 



ipfit) = exp 



1 



cosh t sinh t 



- 1 



exp 



2t 



sinh(2t) 



- 1 



(6) 



Probability distributions corresponding to ip^, ip§ and ipf are infinitely divis- 
ible with finite logarithmic moments. 

Proofs follow from (4) and the properties (c) and (d) of the selfdecom- 
posable distributions. 

Let us note that ipf is the characteristic function of the compound Poisson 
distribution with summand being the sum of independent rv's with the cosh 
and sinh characteristic functions. 

Finally, on the level of the Levy measures N of Y(l), from the BDLP's 
in the property (d), we have the following : 



N d (-) = N § (-) + N f (-), where 

dN § (x) _ 7T 1 



dN<j(x) n cosh( 



7TX N 

2 ' 



dx 



4 sinh 2 )' 



dN f (x) 



TT 



1 



8 cosh 2 



dx 4sinh 2 (f) dx 
Explicitly, as in (4), on the level of the BDLP one has factorization 



(7) 



exp[— ttanht] = exp[l — tcotht] • exp 



t 



cosh t sinh t 



- 1 



(8) 



Taking into account all the above and the property (d) we arrive at the 
identities : 



f , 7r cosh(7ra;/2) 

/ (1 — costx) 2 dx = ttanht; 

Jr\{o} 4 sinh (ttx/2) 

L 



(1 — costrr) 



(1 — coste) 



R\{0} 



4 sinh 2 (7ra;/2) 

TT 1 

8 cosh 2 (7rx/4) 



dx = t cotht — 1; 

2t 



dx 



1 - 



sinh 2t 



(9) 
(10) 
(11) 
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Furthermore, since ipf corresponds to a compound Poisson distribution, the 
last equality implies that 



where we recover the known relation between (cosh w)~ 2 being the probability 
density corresponding to the characteristic function si °£ at and vice versa by 
the inversion formula; cf. P. Levy (1950) or Pitman and Yor (2003), Table 6. 

4. Stochastic interpretation of BDLP's for hyperbolic functions. 

The functions V'cW anc ^ were identified as characteristic functions of 

the background driving random variable Y(l) (in short: BDRV) for cosh and 
sinh SD rv in Jurek (1996), p. 182. [By the way, the question raised there 
has an affirmative answer. More precisely: (23) implies (24). To see that 

note that D 1 = + C, where Di is a copy of D 1 and independent of C . 
The notations here are from the paper in question.] 

More recently, in Jurek (2001) it was noticed that the conditional char- 
acteristic function of the Levy's stochastic area integral is a product of the 
sinh characteristic function and its BDLP ipg. Similar factorization one has 
in Wenocur formula; Wenocur (1986). Cf. also Yor (1992a), p. 19. 

In this section we give some "stochastic" interpretation of the character- 
istic functions ip^(t) and ip§(t) in terms of Bessel processes. 

Let us recall here some basic facts and notations from Pitman and Yor 
(1982) and Yor (1992a, 1997). Also cf. Revuz and Yor (1999). 

For 5- dimensional Brownian motion (B t ,t > 0), starting from a vector a, 
we define the process X t = \B t \ 2 ,t > 0, which in turn defines the probability 
distribution (law) Q s x ,x := \a\ 2 , on the canonical space f2 := C([0, oo); [0, oo)) 
of non-negative functions defined on the half-line [0, oo), equipped with the o~- 
field T such that mappings {u — > X s (u)} are measurable. In fact, (X t , t > 0) 
is the unique strong solution of a stochastic integral equation 



where (/3 t , t > 0) is 1-dimensional Brownian motion. 

The laws Q s x satisfy the following convolution equation due to Shiga- 
Watanabe: 




(12) 




Qt*Qt> = Q[ 



,5+8^ 

x+x 



for all 5, 5 , x, x > 0, 



(13) 
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where, for P and Q two probabilities on (Q, J 7 ), P-kQ denotes the distribution 
of (X t + Y t , t > 0), with (X t , t > 0) and (1^, t > 0) two independent processes, 
respectively P and Q distributed; cf. Revuz and Yor (1999), Chapter XI, 
Theorem 1.2. 

Similarly, let be 5-dimensional squared Bessel bridge of (X s , < 

s < 1), given X 1 = y, viewed as a probability on C([0, 1], [0, oo)). 

Below we use integrals of functionals F with respect to measures Q over 
function spaces. To simplify our notation, as in Revuz and Yor (1999), we 
use Q(F) to denote such integrals. From Yor (1992a, Chapter 2), the Levy's 
stochastic area formula is given in the form 

QLo[exp(-y jf 'dsX.)] =(-^) V2ex p(-^(^ c ^hA-l)) (14) 

However, since Q 5 X ^ = Q S ^ *Q° X ^ , (cf. Yor (1992), Pitman and Yor (1982)) 
we have in fact that 

A 2 x 
g°_ (exp(-y J dsX s )^ = exp ( - | (AcothA - 1)) 

Thus we may conclude the following 

COROLLARY 2. The BDLP Y, for the SD characteristic function <f> d (t) = 
■^^l, is such that Y(l) has the characteristic function 

V>s(0 = exp(l - tcotht) = <5Lo( ex P(^7 (/o i ds x s ))) G J Ao g , (15) 

where (j s , s > 0) is a Brownian motion independent of the Bessel squared 
process X. 

[Here it may be necessary to enlarge the probability space to support in- 
dependent 7 and X .] 

In a similar way, in view of Yor (1992a), Chapter 2, we have 

sf A 2 f 1 \ / 1 \<V 2 x 

Qi( eM- Y I isX,)) = (-j^) exp(--A tanhA), 
so, in particular, 

<5°(exp(-y j dsX s )j = exp ( - ^AtanhA^. 
Thus, as above, we conclude the following: 
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COROLLARY 3. The BDLP Y, for the SD characteristic function ^- t , 
is such that Y(l) has characteristic function 

ip d (t) = exp(-ttanht) = Q° 2 ( exp(it~f {f i dsXs) )) G ID log , (16) 

where a process (j s , s > 0) is a Brownian motion independent of the Bessel 
squared process X . 

Let us return again to functions ip^it) and tp§(t), given in (6), but viewed 
this time as Laplace transforms in t 2 /2. From Yor (1997), p. 132, we have 



1 



cosht 



E 



ex P (-^T«) 



t 



sinht 



E 



exp(-^Tf) 



(17) 



where T{ := inf{t : lZ t = 1} denotes the hitting time of 1 by 5-dimensional 
Bessel process Tif\ t > 0, starting from zero. Jeanblanc-Pitman-Yor (2002), 
Theorem 3, found that the corresponding BDLP's Y are of the form 



Y(h) 



dul (K^<iy 



h>0, 



(18) 



where (r^, h > 0) is the inverse of the local time of TZu^ at r; cf. Revuz 
and Yor (1999), Chapter VI, for all needed notion and definitions. From the 
above we also recover the formulae 



E 



A 2 



exp( — — / dul(\ Bu \<i)) = exp(-AtanhA), 



E 



exp(— 



.0) 



exp(-A(cothA- 1)). (19) 



These as well provide another " stochastic view" of the analytic formulae for 
the BDRV of two SD hyperbolic characteristic functions in (4), i.e., 1/ cosht 
and t/ sinht. 
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